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ABSTRACT 

A detailed study is presented of classical field configurations describing nonabelian 
vortices in two spatial dimensions, when a global SO (3) symmetry is spontaneously 
broken to a discrete group IK isomorphic to the group of integers mod 4. The vortices in 
this model are characterized by the nonabelian fundamental group tti{S0{3) /K) = Qg, 
which is isomorphic to the group of quaternions. We present an ansatz describing 
isolated vortices in this system and prove that it is stable to perturbations. Kinematic 
constraints are derived which imply that at a finite temperature, only two species of 
vortices are stable to decay, due to 'dissociation'. The latter process is the nonabelian 
analogue of the well known instability of charge Iq'I > 1 abelian vortices to dissociation 
into those with charge |g| = 1. The energy of configurations containing at maximum 
two vortex-antivortex pairs, is then computed. When the pairs are all of the same 
type, we find the usual Coulombic interaction energy as in the abelian case. When 
they are different, one finds novel interactions which are a departure from Coulomb 
like behavior. These results allow one to compute the grand canonical partition function 
(GCPF) for thermal pair creation of nonabelian vortices, in the approximation where 
the fugacities for vortices of each type are small. It is found that the vortex fugacities 
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depend on a real continuous parameter a which characterize the degeneracy of the 
vacuum in the model considered. Depending on the relative sizes of these fugacities, 
the vortex gas will be dominated by one of either of the two types mentioned above. In 
these regimes, we expect the standard Kosterlitz-Thouless phase transitions to occur, 
as in systems of abelian vortices in 2-dimensions. Between these two regimes, the 
gas contains pairs of both types, so nonabelian effects will be important. The results 
presented should provide a starting point from which to investigate screening properties 
of the nonabelian vortex gas. 
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1 Introduction 



The study of vortex defects in both two and three spatial dimensions, has revealed many 
interesting properties which have had a wide variety of applications in both condensed 
matter physics [1], and cosmology [2]. In their guise as cosmic strings [1], vortices in 
three spatial dimensions, which occur as the result of the spontaneous breakdown of 
a global symmetry in the early universe, may go some way in explaining the origin 
of large scale structure. At the other extreme, it has been known for some time [3] 
that vortices in 3-dimensional condensed systems such as superfiuid helium IV, play 
an important role in understanding the properties of these systems. 

It is well known that vortex defects in any dimension arc characterized by 7ri(M), 
the fundamental group of the vacuum manifold M of the theory. The type of vortex 
defects that occur in for example helium IV are 'abelian' in the sense that they are 
characterized (as we shall review later), by an abelian fundamental group, isomorphic 
to the group of integers Z. In 2-spatial dimensions, abelian vortices have particularly 
simple interactions which allows one to write down the grand canonical partition func- 
tion (GCPF) for thermal pair creation to all orders in the vortex fugacity in terms of the 
so called Coulomb gas model [4] . The pioneering work of Berezinski [5] , and Kosterhtz 
and Thouless [6], into their statistical mechanical properties, showed that a gas of such 
vortices underwent a novel kind of phase transition at some critical temperature Tp. 
A simple physical picture emerged where for T < Tc, vortices and antivortices form a 
medium of bound pairs which subsequentally dissociates into free vortices and antivor- 
tices for T > Tc- These results when apphed to approximately 2-dimensional systems 
such as helium IV thin films, lead Kosterlitz and Nelson [7] to predict a universal jump 
in the superfiuid density at Tc, which was later experimentally verified [8]. 

In this paper we shall study some properties of nonabelian vortices in two spatial 
dimensions, which should eventually lead to an understanding of their statistical me- 
chanical properties. There is particular interest in these kind of vortex systems, not 
only because they are interesting extensions of those described above, but also because 
they are known to occur in certain 2-dimensional liquid crystals [9]. In fact, we shall 
consider the simplest model of spontaneous global symmetry breaking that yields a 
nonabelian 7ri(Af) = Qg where the latter is the order 8 group associated with the field 
of quaternions. Vortices described by this fundamental group are just those which 
are expected to occur in so called nematic liquid crystals [9]. As we shall see in later 
sections, the interactions of nonabelian vortices arc very much more complicated than 
their abelian counterparts, so much so that we cannot expect to be able to write down 
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the complete nonabelian vortex gas partition function to all orders in the relevant fu- 
gacities. We will only calculate the latter to fourth order in fugacities. Whilst this 
may seem very far from the all orders results in the abelian case, it is still sufficient 
to investigate whether a kind of nonabelian charge screening occurs, as we shall now 
explain. 

Ordinary abelian charge screening is the mechanism behind the Kosterlitz-Thouless 
type phase transition described above [1]. There are a number of ways of exhibiting this 
screening, for example by exploiting the map between the vortex gas partition function 
and that of sine-Gordon theory in 2-dimensions [10]. Alternatively one may derive 
a so called Poisson-Boltzmann (P-B) equation [11], satisfied by the linearly screened 
potential between a test vortex and antivortex placed in the gas. The latter equation 
is derived in a perturbative small fugacity expansion, and it is clear that in order to 
compute the first order screening effects one needs to know the energy of a configuration 
with at least two vortices and two antivortices. One of these vortex-antivortex pairs 
will play the role of test charges, whilst the other will describe lowest order screening 
due to thermal pair creation. Of course due to the abelian properties, the resulting 
P-B equation is relatively simple as are its solutions [1]. One learns however, that it 
is not necessary to have an exact expression for the partition function, to investigate 
screening properties. 

From the results given in this paper, it should be possible to generalize the no- 
tion of linearly screened potential, and to derive a nonabelian version of the Poisson- 
Boltzmann equation [12]. 

The structure of the paper is as follows. In Section 2 we introduce the 2-dimensional 
field theory model which describes nonabelian vortices associated with 7ri(M) = Qg. 
After a brief review of the homotopic classification of vortices, we present an ansatz 
describing isolated nonabelian vortices, and give numerical minimum energy solutions 
for the Higgs fields that define their 'core' region. The properties of these solutions are 
shown to guarantee finite energy. In Section 3 we present an ansatz for multivortex 
configurations and write down an expression for the GCPF describing the thermal cre- 
ation of vortex-antivortex pairs. The energy density of all relevant configurations up to 
a maximum of four vortices is calculated. Section 4 describes a method for obtaining 
the exphcit form of the total energy of each of these configurations. As an illustrative 
example, the method is used to compute the energy of a particular configuration in 
which nonabelian effects are expected to be important. After a conclusion and com- 
ment on the results, Appendix A gives details of a proof that our ansatz for isolated 
nonabelian vortices is indeed one of minimum energy. We should say at this point that 
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there have been a number of recent papers concerning the subject of nonabehan vortices 
(e.g. ref.[13]). However these differ from our approach in that the symmetry breaking 
leading to nonabehan vortices is local. Moreover, the main investigations of these pa- 
pers concerns the quantum properties of nonabehan vortices and their connections to 
systems obeying fractional spin statistics. 

2 Definition of the model 

In this section we shall describe in detail, a Higgs system which will give rise to non- 
abehan vortices based on the discrete group Qg of quaternions. As stated previously, 
the latter are probably the simplest examples of nonabehan vortices, but which also 
have direct physical interest. We shall consider the Higgs model first described in [14], 
where the scalar field or order parameter $ is in the five dimensional representation 
of 5*0(3) , and so can be written as a traceless, symmetric 3x3 matrix. The energy 
density of the model is taken to be 



and gab is the 2-dimensional flat euclidcan metric, where a, b = 1, 2 run over 2- 
dimensional polar coordinates (r, 6 ) respectively. 

Any field configuration $ can be diagonalized by some orthogonal transformation 



S = -Tr [g'^' (9,$) {dt,^)\ - V{^) 



(2.1) 



where the potential is chosen to be 




(2.2) 




(2.4) 



(2.3) 



The extrema of the potential, are given by the solution to the equations 
= A(2(^^ + ^2 + '^^1^2 + 3</7i</?2) + P{—vl — 2¥'l</?2) 



(2.5) 
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- + 2^2) + \'{2^\ + A^l + Q^l^2 + Q^i^l) 

(2.6) 

Writing ip2 = d'-Pi for some real parameter a, we have to solve two simultaneous cubic 
equations ( p.5|) , ( p.6|) in say the variable y^i. For p = 0, but arbitrary A, A' one finds 
that there are degenerate set of minima which lie on an ellipse in the ^pi — ip2 plane, 
defined by the equation 

+ ^^^^+ ^l) - 1] = (2.7) 

where in eqn.( p.7|) , and in the rest of the paper, we absorb the mass scale /iq into the 
fields ipi, ip2. For varying values of A and A', eqn.( |2.7| ) defines a continuum set of 
ellipses, illustrated in Fig.l . If p 7^ 0, then in fact this degeneracy is broken to a large 
extent, and the potential has three absolute minima at the points ^2 = o-Vi the 
same ellipses, with a = 1, —1/2 and —2 respectively. These minima are indicated by 
the intersection of the three dashed radial lines with the ellipses in Fig.l. Since the 
vacuum expectation value < $ > can be written in terms of (pi and a at a minimum, 
one has 

/I \ 

<$>=(/?! a (2.8) 
-l-aj 

It is clear from eqn. (|2.8| ) that when p 7^ 0, the three absolute minima correspond to an 
unbroken f/(l) global symmetry, when < $ > is invariant under rotations about the 
X, y and 2;-axes , for a = 1, —1/2 and —2 respectively. When p = there is again an 
unbroken U{1) at the above three points as well as at three new points corresponding 
to the intersection of the ellipse with the dotted radial lines in Fig.l. The dotted 
lines shown in this figure are just the continuation through the origin of the dashed 
radial lines, which together define the lines = O'^i for the three values of a just 
given. For p 7^ 0, the vacuum manifold M, i.e., the set of points in the space of 
matrices given by $ = ^ < $ > Q^^ , Q G 5*0(3) , will be given by the quotient space 
M = S0{3)/U{1). The fundamental group of M in this case will be isomorphic to 
^2, the group of integers modulo 2. As such, the vortices which can arise after this 
symmetry breaking will necessarily be abelian in nature. Since we wish to investigate 
vortices of the nonabelian variety in this paper, we shall take the coupling p always to 
be zero. Then as long as we choose values for < > and < (p2 > on the ellipse which 
do not coincide with the six points at which U{1) symmetry remains unbroken, the 



5 



maximum symmetry group of < $ > will be a discrete subgroup IK of 5*0 (3), which is 
isomorphic to the group of integers mod 4. As we shall see in more detail below, this 
leads to the manifold M having a fundamental group given by the nonabelian discrete 
group Qg, isomorphic to the group of quaternions. Throughout the rest of this paper, 
we shall concentrate on this case. 

We now move on to discuss in detail an ansatz for isolated nonabelian vortices with 
specific reference to the group Qg. In appendix A we shall show that our ansatz is the 
correct one in that it describes vortices with the lowest energy. The vacuum manifold 
for the case when < $ > breaks 50 (3) down to the discrete subgroup IK which has 
elements 
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1/ 



is given by M = 5'0(3)/lK = S'f/(2)/Qg, where we have used the well known ho- 
momorphism between 5*0(3) and SU{2). It follows that vortices are characterized by 
7Ti{M) = Qg. Using the homomorphism between SO (3) and SU{2) we can equally 
work with fields $ which are representations of either group. We will usually work 
with 5*0(3) representations and group elements (the latter denoted by Q ) but will 
often give the corresponding SU{2) fields and elements (denoted by (7 ), on which Qg 
has a simple action. 

Our ansatz for a single vortex at the origin in the (r, 6) plane will be 



(2.10) 



where in eqn.( |2.lOD , ^^diag denotes the diagonal matrix (v5i(r), (p2{''^)i ~(v^i(^) + V'2('")) )• 
It is assumed that v'i,v52 are functions of r only whilst the group elements Q depend 
only on 9. It is important that on encircling the vortex at r = 0, i.e. as 6 is varied 
from to 277, $(r, 9) be strictly single-valued. The corresponding boundary conditions 
on G{0) consistent with this constraint are. 



g{e + 2Ti) = g{e)h (2.11) 

where h is an element of K. It might appear that there are only four 'types' of vortices 
corresponding to the four elements of K. However it should be remembered that 50(3) 
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is not simply connected since ni{S0{3)) = Z2, which effectively doubles the number 
of homotopically inequivalent closed curves on M. This is more transparent when 
working with SU{2) group elements. Under the homomorphism SU{2) = 5*0(3), the 
group of quaternions Qg gets mapped to the abelian group K. So somewhat hidden 
in the boundary conditions given in eqn. (|2.11|) is the fact that closed curves on M are 
actually characterized by elements of Qg. 

Next we need find out what are the explicit representations of 7ii{M). From a 
mathematical point of view there are of course, an infinite number of homotopically 
equivalent representations. However, we are only interested in the physically relevant 
ones, i.e. those which yield a minimum energy for the corresponding vortex, and we 
expect that such 'minimal' representations will be unique. The elements K,2-, ■■■■,K,i 
given in eqn.(|2.9|) correspond to discrete rotations by 180 degrees about the x^y and 
2;-axes respectively. Therefore, a simple ansatz for the S'0(3) elements Q which satisfy 
the boundary conditions in eqn. (|2.11| ) for a given element of K, is for them to be 
rotations by angles 9/2 about the respective axes. In fact, we shall show in Appendix 
A that this ansatz actually corresponds to minimal energy vortices. 

In Table 1, we have listed these elements, together with the corresponding SU{2) 
representations. Also listed are the discrete group elements that appear in the boundary 
conditions as we encircle the vortex. The SU{2) representations in particular, make 
it clear that this group is Qg. We have used the notation Qi, Qj, Qk, G^^, Qj^ , Qk^ and 
to denote group elements whose boundary conditions are twisted by the elements 
j, /c, — i, — j, —k and -1 of Qg, where 

^'^=f = k'' = -\ 
ij = k , jk = i , ki = j (2.12) 

are the multiplication rules of Qg. The parameter 6q in Table 1 is an arbitrary constant 
associated with the origin of the 6 coordinate. 

Given these explicit forms for Q, we want to minimize the energy in eqn.( ^.l| ) in 
the remaining fields (pi{r) , ip2{r) , in order to find stable vortex solutions. Because of 
the $ —>■ — $ symmetry of the energy, (assuming p = ), the equations for minimizing 
the energy will be the same for Qa as for G^^, for a e Qg. We shall refer to these as 
being vortices and antivortices of type a respectively. Furthermore, the vortices of type 
-1 are by nature abelian since -1 commutes with all other elements in Qg. Moreover, 
we shall see in the next section that in a finite temperature system they are unstable 
to dissociation into either i,j or k type vortex-antivortex pairs. For these reasons we 
will concentrate for the remainder of the paper, on vortices of type i,j and k. The 
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parameter a which denotes species type, will therefore implicitly run over only these 
particular types. 

Substituting the ansatz for Q^^ in these three cases into the energy we get 



2 



- V 



JO 



+ J-P^{^,,^^)-V{^i,V2)] (2.13) 



2r 



The effect of the terms Qa{0) on the energy is to add certain potential terms, 
denoted by Paifi, V2) in eqn.( |2.13| ) which are given by 



Pi = + 2 ¥^2)' (2.14) 

P, = ^(2^1 + ^2? (2.15) 

Pk = ^(^1 - V2? (2.16) 
The minimization equations for Ea are 

r r ir'^ OLpi Oifi 

2a„V92 + - dr(p2 + d„Vi + -dripi - TT^^f^ + = (2.17) 

r r / r"' 0(p2 c'V^2 

On physical grounds, (namely the finiteness of the core energy E^, which we will be 
defined later), we anticipate that the solutions to these equations that minimize the 
energy have the following asymptotic behavior in the radial coordinate r. As r — 0, 
it is sufficient that (fi and (f2 both vanish at least as fast as r, whilst for r oo, we 
expect both (pi and ip2 go to their vacuum values, which minimize V. To show that 
smooth interpolating functions exist with these asymptotic properties, one may solve 
the non-linear differential equations ( |2.17| ) numerically. This was done by discretizing 
the energy in eqn.( |2.1|) and carrying out numerical minimization, for various values of 
the parameters A, A' and vacuum parameter a. Fig. 2. illustrates the solutions found in 
the case of an z-type vortex with A = 0.01 and 2, A' = and a = —1 (we shall explain 
why we are interested in the value a = -1 later). In Fig.2, the dotted curves show the 
numerical solutions whilst the solid ones are those of a trial solution, needed as input 
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in the minimization procedure. These trial solutions are of the form 

1 f rfxo \ 



A(l + a + a2) V'^/^o + l, 



trial 



^2\tr.al = , (2.18) 



The normalization factors in eqns. (|2.18|) are such that (y9i(oo), ^92(00) lie on the ellipse 
of degenerate vacua, described in the previous section. 

As one can verify from Fig. 2, these trial solutions are very good approximations, 
particularly close to the origin of r. It is also apparent that there is a 'core' region 
(marked as r = in Fig. 2) beyond which the values of y^i, do not change appreciably, 
from their vacuum values 0. Of course, different definitions of are possible, depending 
on how close to unity one requires the ratios ipi/ < ipi> and ip2l < ^2 > be. The 
lack of a precise core radius is not a particular property of nonabelian vortices, but is 
also apparent in the abelian case. What is clear in both cases however, is that different 
choices for Tc will simply give rise to constant factors in the expression for the self 
energy of a vortex (see eqn.( |2.l9| ) below). Similar solutions to these exist for j and k 



type vortices, the only difference is in the algebraic form of the potentials Pa which, 
it turns out, has a minimal effect on the behaviour of the solutions. It will be a very 
useful (and good) approximation to replace the numerical solutions for ^\,^2 of ^5 3 
and k type vortices by those of the trial solutions given in eqns. ( 2^181 ) for the region 
r < re-, and by < if)\> ,< ^2> for > ^c- 

This being the case, it is straight forward to compute the energy of a single i, j or 
k type vortex (or antivortex), placed in a 2-dimensional circular box of radius R. The 
result is 

Ea = E^ + n <Pa>ln (—) a = k (2.19) 



•.re 

where is the core radius of the corresponding vortex configuration, and is typically of 
order unity in units of inverse mass /ig ^. The logarithmic self energy terms in eqn. (|2.19|) , 
correspond to contributions coming from the region r > Vc- The core energy, E^ in 
eqn.( |2.19D can be approximated by taking the range of integration in eqn.( |2.13D to be 

^ Although there is no definite value to rc, we can define it to be that value of r such that 
ifi/ < (fi >= (^2/ < ^2 >= 1 ^ ■^i that is ifii, f2 are within ^ of their vacuum values. This gives 
flora ^ 1-72 
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between r = and r = Tc, and is given by 

rrc / r \ ^ 

= vr <P, >{ln(r, + l) + ^--l} + /(<(^i>,<y.2>,re) (2.20) 



The first term in eqn.( p.20D comes from tlie potential Pa present in due to vortices 
of type a. Tlie second term /(< ipi >,< >,rc) represents contributions coming 
from other terms in ( |2.13D , which are independent of the vortex type a. These can 
be easily calculated, but exphcit expressions will not be needed in what follows. It is 
worth comparing the results of eqns.f pTT^ ), (|2.20|) to those one would obtain for a single 
abelian vortex in 2-dimensions, such as those produced in spontaneous breaking of U (1) 
symmetry [15]. In the abelian case we would also obtain an energy formula similar to 
eqn.( |2.19D , with the coefficient of the logarithm term proportional to < 0^ >, where 
q is the vorticity and < > is the U{1) breaking vacuum expectation value. Hence the 
factors of < Pa > in eqn.( |2.19| ), ( p. 201 ) for fixed < ipi > and < ip2 > are proportional 
to a kind of nonabelian vorticity or nonabelian charge. The difference between the 
abelian and nonabelian example is that in the former, we can have arbitrarily large 
(but quantized) vorticity q, whilst in the latter there are clearly only a finite number 
of possibilities, given by the range of a. This is a consequence of vri(M) = Z, ( Z 
being the group of integers), and vri(M) = Qg in the abelian and nonabelian cases 
respectively. 

Finally, we can understand the logarithmic behavior in eqn. (|2.20|) and in particular 
why there is a logarithmic divergence as R —>■ oo, even for nonabelian vortices. This is 
because, as the numerical results suggest, one can think of these vortices as having a 
definite core region, which for distances r 3> Vc, appears approximately as a point-like 
charge. It is natural then that the field energy outside such an isolated point charge 
be divergent as we go to the infinite volume limit, whilst the logarithmic behaviour 
is simply a consequence of the 2- dimensionality of the system. In fact, to see the 
truly novel effects of nonabelian vortices we have to study more than simply self- 
energies. Clearly we have to consider multi- vortex configurations where at least two of 
the vortices are of types h and g where the latter are two non-commuting elements of 
Qg. We shall study this situation in the next section . 
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3 Systems of multi-vortices 



Having discussed at length in the previous section the ansatz corresponding to isolated 
nonabelian vortices we now move on to discuss multiple vortex configurations. The 
simplest such configuration would be a vortex-antivortex pair of type a. Since we 
are ultimately interested in the statistical mechanical properties of such vortices, one 
should consider what happens in a real system at finite temperature. The logarithmi- 
cally divergent self energy of a single vortex discussed in the last section, prohibits, 
in the thermodynamic limit, single vortex creation in these systems due to thermal 
fiuctuations. This is a property which is familiar in the case of abelian vortices [1]. 
However, neutral vortex-antivortex pair creation is not suppressed in this way, and one 
should therefore expect these to be present at finite temperature. 

Intuitively we expect that whatever the correct description of such a pair is, it 
should go over to the previous ansatzes for an isolated vortex and antivortex in the 
limit where we take the distance between the pair very large with respect to their 
typical core size Tc- As such we shall take as our ansatz for a vortex of type a at the 
point z\ and an antivortex of the same type at Z2 to be 

Q^{Z, Zi, Z2) ~ Ga{z, Zi)g_aiz, Z2) (3.1) 

where in eqn.( |3.1|) , the Qa correspond to the 1-vortex solutions obtained in the previous 
section. The corresponding field configuration denoted by $q, is 

= Ga{z, Zi) g~^{z, Z2)^diag{z, Zi, Z2) Ga{z, Z2) G~^{Z, Zi) (3.2) 

Again, for well separated pairs we expect that the 'core' field ^diag vanishes as z 
approaches either of the vortex centre's at zi or Z2 , (which ensures a finite core energy 
in each case) just as in the 1-vortex case. In fact this behavior can be confirmed if one 
were to solve the minimum energy equations for ipi, ip2 for the vortex pair with given 
by eqn.( |3.2|) . As z approaches either of zi or Z2, the form of Qa goes over to that of the 
single vortex, modulo an unimportant constant shift in the angular coordinate 6. Hence 
from our results of the last section, it is clear that the solutions for (pi, ip2 will become, 
in the limit of large pair separation, those described in the previous section. Thus 
as before we can approximate ^pi.,^2 to be given by their vacuum values at distances 
\z — Zs\ > Tc, where s = 1, 2 . For distances \z — Zs\ < we can approximate the 
fields by the trial solutions given in eqn.(^7TB|). Putting all this together we find that 



the energy of a vortex-antivortex pair is given by 

E^{z^, Z2) = 2EI + 2El^'f + E^*(zi, Z2) (3.3) 
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where in eqn.( p.3| ) the core energy E^^ was given in ( |2.20| ). The self energy E^^^ and 
interaction energy i?^"* are computed to be 

i?-^/ = vr < > In(-) 

E^\z,,Z2) = 27r <P„>ln|^i-^| (3.4) 

K 

where R is the size of the system, and the functions Pa were defined in eqns.( PTT^ )- 
(2.17). Note that we can take the thermodynamic hmit i? — > oo in eqn.(2]4| ) since 
R cancels between the self and interaction energies. This is simply a consequence of 
taking a vortex -antivortex pair-which can be thought of as neutral^. This phenomenon 
is already known to occur in the case of abelian vortices, and indeed E^^^{zi, Z2) above, 
is just the standard Coulomb type of interaction which one also finds in this case. This 
is not surprising because simply studying a vortex-antivortex pair of type a does not 
probe the nonabelian structure of the theory since Qa and Q-a are commuting elements. 
Clearly to study the effects due to non-commuting elements Q^^Qis, we have to look 
at least to configurations containing four vortices: a vortex-antivortex pair of type a 
interacting with a pair of different type /3. 

Before studying the details of this, we can make life easier by consideration of what 
one might call kinematical constraints associated with the thermal creation of vortices. 
As we shall see this implies that for any given value of the vacuum parameter a, only 
two types of vortices will be stable to dissociation into one another, where these types 
depend on the specific value of a. 

Let us imagine creating vortex-antivortex pairs in a system by gradually increasing 
the temperature. The self energy eqn.( p.4|) of a vortex of type a in the pair is pro- 
portional to the quantity < Pa > (see eqn.(p.4|). As we discussed earlier, the < Pa > 
factors can be thought of as being proportional to a generalization of vortex charge or 
vorticity to the nonabelian case. As such one has to consider whether a nonabelian 
vortex of given charge is unstable to 'splitting' into two other vortices in such a way 
as to conserve overall charge, yet lower the self energy. 

To understand this process let us consider the abelian case once more. Here the 
vortex self energies are proportional to , where q is the quantized vortex charge. 
Thus it is clear that a charge q = ±2 vortex will be unstable to dissociation into two 
charge ±1 vortices because -^^=±2 > 2Pgl'£]^. The same is true for all higher charged 
abelian vortices with |g| > 1, so that in this sense If?] = 1 vortices are ultimately stable 

In the nonabelian case, this means that G{z, zi, Z2)is single valued upon simultaneously encircling 
the points zi and 22- 
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[15]. This explains why only q = I vortices are considered in the context of Coulomb 
gas models and K-T type phase transitions [1]. We want to apply a similar analysis 
in the nonabelian case. The main differences are that we have a finite number of 
possible charges in this explained earlier, and the rules for possible dissociation 

of vortices are governed by the group relations of Qg not simply by the addition of 
charges. 

It should be stressed at this point that the stability to dissociation we are discussing 
is different to, and must be considered in addition to, the perturbative stability analysis 
performed in Appendix A. Clearly the process of dissociation of a vortex is non-local, 
and is separate from the stability of a vortex to local perturbations. 

To begin with consider the self energy of a -1 type vortex. We saw in the previous 
section that there were three kinds of ansatz for the 5*0(3) element, denoted by ^i^i, 
where ( runs over the elements i,j,k of Qg (see Table 1). It is easy to see that the 
corresponding self energies E^l'^'^'^^ satisfy the relations 

^-^/,(C) ^ 4^-«/ > 2Ef^ (3.5) 

for arbitrary values of a. It is clear that -1 type vortices are always unstable to 
dissociation into a pair of vortices of type i,j or k. Which actual decay branch occurs 
depends on the specific value of a and will be whatever is the lowest energy of the three 
possibilities E'^_^l^^'^'^\ From a practical point of view however, this result means that 
we need not consider -1 type vortices at all in investigating the statistical mechanics of 
nonabelian vortices. They are like the Iq'I > 1 abelian vortices referred to above. Let us 
now consider the remaining cases of type a (clearly since Eff = E'_'!/ the arguments 
will hold equally for antivortices ). 

Using the group composition rules for Qg, a A; type vortex for example, can dissociate 
into an i and j type (or -i and -j type ) if 

{1-af > {l + 2af + {a + 2f i.e. - 2 < a < (3.6) 

where in eqn. (|3.6| ), we have substituted for in terms of the parameter a. Thus 
in the interval —2 < a < — 1/2 , referred to as region lij in Fig.3 , k -type vortices 
dissociate, and at most we would expect to just find vortices and antivortices of type i 
and j in this region. At the same time we learn that i and j -type are stable to further 
dissociation, since in either case the path (restricted by the composition rules of Qg), 
necessarily involves producing a k type vortex. Notice that while lij is a connected 
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portion of the real line of a values, the corresponding regions in the space of vacuum 
expectation values < (fi >, < (f2 > are two disconnected pieces. A similar analysis of 
possible decays of i and j types for a values outside of region lij gives the following 
results, 

Eff > Eff + Eff 

{l + 2af > {2 + af + {l-af i.e. a < -2 or a > 1 (3.7) 



and finally 



j k i 



+ > {l + 2af + {l-a)^ i.e. - - < a < 1 (3. 



1 
2 

the different regions of a given in eqns.(|3^, |378| ) are denoted by Ijk and Iki in Fig.3. 



Fig. 3 clearly shows that the three regions fit together exactly, to cover the entire ellipse 
of vacuum states, without overlap. This means that the two types of vortices indicated 
by the subscripts on and I^i are stable to further dissociation. It is interesting 

to see that all regions are separated by points with either a = —1/2, 1 or 2, which 
correspond to an unbroken U{1) symmetry (see Section 2 ). At such points there are 
no vortices, so it appears that one cannot smoothly go from one region to another by 
varying a. In fact the relative positions of the regions Iij,Ijfc and 1^, in Fig.l. can be 
understood in the following way. 

The point where any two regions touch has an unbroken U{1) symmetry that cor- 
respond to rotations about either the x,y or z axis, depending on which of the six 
points, where regions lij etc meet, we consider. For example at the point where lij and 
Ijk meet, the U{1) symmetry corresponds to rotations about the y-axis (which one can 
easily see from the form of < $ > (see eqn.( p.8| )) at this point). In some sense, we can 
now understand why this particular symmetry should appear here because in 'passing' 
from lij to Ijk, the vortices of type i have been rotated by 180 degrees about the y-axes 
into those of type k whilst leaving those of type j invariant. The other regions in Fig.3, 
and their relation to the particular U{1) symmetry present at their intercept, can be 
similarly understood. 

The end result of the above kinematic considerations, is that for any given ra- 
tio < <y92 >/ < fi >= CL, we need only concern ourselves with the statistical mechan- 
ics of two species of vortices. Let us for definiteness choose a to be in the region 
— 2 < a < —1/2, i.e. we only study i and j type vortices. 

Within this approximation, we can write down a general expression for the grand 
canonical partition function (GCPF) Z of thermal pair creation of the above species 
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of nonabelian vortices, 



Ni=ONi=0 Ui=lu,=l-^ "'c "'c 



X 



E E exp|--E(iV„iV,,6;.,,a,,^.„^.,)} (3.9) 
{e.J {6.,} 

n 5,^=1 n 5,^=1 

In eqn.(|0| ), the fugacities Zj,Zj associated with vortices of type i and j, are defined 
in terms of the core energies E'^ given in eqn. (|2.20| ) 



z 



'3 



exp(-E,7T)) 

exp(-E|/T) (3.10) 



respectively. Q^^ and ^^.^ denote the 5*0(3) group elements corresponding to all the i 
and j type vortices of a particular configuration having total energy E. As we have 
already seen in this section, only configurations with H^Li Qu^ = 1 and Quj = 1 
will contribute to the GCPF above, in the thermodynamic limit. Although eqn. (|3.9|) 
is a complete formula for the GCPF, we want to study it perturbatively in a small Zj 
and Zj expansion since we do not know how to compute E for arbitrary configurations 
of i and j type vortices. 

Having considered in some detail the case of single vortex -antivortex pairs, we 
now move on to consider the next order terms (in a double expansion of Z in the 
fugacities Zj, Zj ). These will be configurations containing four vortices and antivortices 
arranged in two pairs, because the overall neutrality condition mentioned previously, 
demands that we consider only even numbers of vortices and antivortices. The simplest 
configuration of this type consists of two vortex-antivortex pairs of the same type i or 
j. Because the pairs are of the same type, they commute with each other, so we expect 
that the resulting energy be that of the corresponding Coulomb system of charges 
[1,4]. The ansatz for this system will be just the generalization of the single pair 
ansatz given in eqn.( |3.1| ), which again we expect to be a reasonable approximation at 
inter-pair distances greater than the core size r^.- Thus with the understanding that the 
index a now runs over i and j types only, the S0{?>) group element, Qa,a, corresponding 
to vortices of type a at points Zi and 2:3 and antivortices of the same type at ^2, ^4 will 
be 

Qa,a{zi, ■ ■ ■ , ^4) ~ Ga{z, Zi) Ga{z, Z2) Qa{z, Z3) ^^(z, 2:4) (3.11) 
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with the corresponding form of $, which we denote by ^a,a being 

y<^diag{z, Zi,Z2,Z3,Z4^) 
XGa{z,Z4)g~^{z,Z3)Qa{z,Z2)g~'^{z,Zi) ] a = i, j 

(3.12) 

Again for the purposes of computing the energy of the configuration given in 
eqn.( p.l2"D , the core field ^diag{z, zi, Z2, z^, Z4) will be approximated by the usual trial 



functions of eqn. (|2.18|) whenever \z — zt\ < Vc for t = 1, 4 or by < ^diag > when the 



argument z lies outside this region. The energy of this configuration E^^a is given by 

4 \z - z\ 

EaAzi,Z2,Z3,z^) = AE'^ + u < P^> J2 IsqtlnC— -) (3.13) 

s,t=l 



where in eqn.( p.l3"l) the charges are given by gi = gs = 1,^2 = ^4 = —1- We 



see that the energy q, of two pairs of the same type a has the usual Coulombic 
form. Moreover, in the GCPF eqn. (|3.9| ), we have to sum over similar configurations 
to those above but with different charge assignments Qs at the points Zg- However, 
the relative energy of these different configurations is trivially related to that given in 
eqn. ( |3.13| ) since it simply involves a permutation of the charge assignments. This is a 
manifestation of the property that all the 5*0(3) group elements in eqn.( p.l^ ) commute. 

At this point it is worth mentioning the fact that whilst one cannot write down 
the energy E appearing in the GCPF Z for arbitrary configurations of i and j type 
vortices, it is possible for those configurations given purely in terms of i or j types 
only. The resulting energy E is just the generalization of eqn.( p.l5| ), corresponding to 
two Coulomb gas systems of i and j type vortices respectively. Viewed in this way, the 
GCPF Z describes how these two different, 'non-commuting' Coulomb gases interact 
with each other. The relative strength of each Coulomb system is controlled by the 
fugacities Zj and zj. One can imagine two extreme cases where Zj ^ Zj or zj ^ Zj, so 
that Z is dominated by i or j types. Each of these situations correspond to a kind 
of abelian limit of Z where we would the usual K-T picture of a bound medium or a 
plasma of vortices of either type. 

Since Zj and zj are expressed in terms of the core energies Ef and Ej (see eqn. ( p. 10] ), 
one has to check that these two abelian regimes can be actually realized. Because we 
have used the kinematic arguments given earlier to focus our attention on i and j types 
only, the vacuum parameter a is implicitly taken to be in the range (—2, —1/2). It is 
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then clear from the a dependence of the core energies (eqn.( |3l6[ )), that for sufficiently 
small temperature T, the abelian limits Zj 3> Zj and Zj ^ Zj correspond to the values 
a —>■ —2 and a —>■ —1/2. For a = -1 in the mid range of (—2,-1/2), Zj = zj and 
we are in the extreme nonabelian regime where both vortex types are equally likely 
to be found in the system. Clearly this region is the most relevant to investigating 
nonabelian effects, and clearly explains why we chose a = -1 in the numerical plots 
given in Fig.2 and Fig. 5. 

Let us continue to study the 2-pair contributions to Z by considering the case when 
the two vortex-antivortex pairs are of different types i.e. one of type i, the other j. We 
expect a very different answer in this case since now the full nonabelian structure of 
the theory should be revealed via the non-commutativity of the group elements Qi, Qj. 
In the same manner as above, the form of $ for a pair of type i at the points zi,Z2, 
and pair of type j at 23, ^4, denoted by will be 



*i j (2:, ^1, ■ ■ ■ , Z4) = Giiz, zi) ^{z, Z2) Qj{z, Z3) ^{z, Z4) 

>^^diag{z, Zi, Z2, Z3, Z4) 
xGjiz, Zi) g~^{z, Z3) Qi{z, Z2) Qi^{z, Zi) 



(3.14) 



In the GCPF we also must sum over other configurations related to ^^Ij by permuta- 
tions of the various group elements, keeping the ordering of the points 2:1, 2:4 fixed. 
However, since Qi and Qj do not commute, the energy of these different configurations 
are not in general, simply related to each other. Also they are much more complicated 
than the simple Coulomb potential we saw previously, which is to be expected since 
the latter is really a property of abelian vortex systems. Because of this complexity it 
is not so straight forward to compute explicitly the energy of each configuration. The 
difficulty evident from the form of the energy density, which we will now discuss. 

For the configuration given by (^fj of eqn. (|3.14| ) the energy density, denoted by 
£i(^_i)j(^_j) , is found to be 

8i(^^iyj(_j) = 2{ifi + 2LP2Y - 6v9i (2 ipi + ip2)sm'^[A{z - Z3) - A{z - Z4)] 
X [A,{z - zi) - A,{z - Z2)] [A,{z - zi) - A,{z - Z2)] 
+ 2{2<fi + if2)^ 

X [A^{z - Z3) - A-Xz - za)\ [As{z - Z3) - As{z - Z4:)] 

(3.15) 

where in eqn.( t3.15"|) , A{z — Zg) = Arg{z — Zs) ,s = 1,...,4, and the subscript z{z) 
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on A denotes differentiation by ^ (^) respectively. In computing Si(^-i)j(^-j), we liave 
suppressed terms involving derivatives of and as well as those involving the 
potential V. They will only contribute to the core energies as we discussed earlier 
(see eqn. (|2.19|) ). It is interesting to observe in this equation, that an effective interaction 
is induced between the pair of type i and j given by the (sine)^ term. As a check on 
the correctness of eqn. (|3.15|) , one may readily verify that goes over to the 

energy density of a single i (j ) type pair if one takes the limit Z4 {zi — > Z2 ) 

respectively. 

Now in principle we would have to calculate the energy of another 255 possible 
orientations of the matrices Qi, Qj and their inverses which corresponds to all possible 
permutations of the 'charge' assignments at the points zi, ...,Z4^. However as we shall 
show below, one can employ various similarity transformations acting on the elements 
Qa to connect the energy of one configuration to another. In this way we will see that 
one need only calculate explicitly the energy density of two other orientations other 
than the one described by ^i^j ■ 

Denote by E[Q, < ^diag >] the energy of an arbitrary orientation of the four matrices 
Qi, Q-i) Qj, Q-j at the points 2:1, . Then for any constant matrix U , we have 

E[ugu-\ < >] = E[g, u < > u-'] (3.16) 

Now consider the particular U matrices listed in Table 2. They generate the transfor- 
mations {—i, —j), and on the i and j types of vortices, and hence 
generate permutations of the QaS. Using these transformations, it turns out that we 
can relate the energy of an arbitrary orientation of QaS to those of just three particular 
orientations; $• j'' (already given in eqn. (CT) ) together with $5^, $5^^ defined by 

^\^^{z,Zi,- ■ ■ ,Z4) = gj{z,Zi)Q-i{z,Z2)Qi{z,Z3)g^j{z,Z4) 

X ^diag{z,Zi,Z2,Z3,Z4) 
xGjiz, Zi) Q-iiz, Zs) Qi{z, Z2) G-j{z, Zi) 

and 

^f^{z,Zi, - ■ ■ ,Zi) = gi{z,Zi)gj{z,Z2)G-i{z,Z3)g^j{z,Z4) 

X ^diag{z,Zi,Z2,Z3,Z4) 
xGjiz, Zi) gi{z, Z3) G-jiz, Z2) Q-iiz, Zi) 

(3.17) 

The energy densities Sij(-j)(^-i) and Sij(_^i)(-j) of these last configurations, (again in the 
region where |-2 — > Tc so that one may approximate ^diag by < ^diag >), are found 
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to be 



= {2 (2 + ip2f + sm^[A{z - Z2) - A{z - Z3)] 
X ( 6 (2 v?i + 992) ^2 sin^ [A{z - Z2)] 

-6</?i {(fii + 2(P2))]Az{z - zi)A^{z - zi) 
-{4 (2ipi + ip2f cos[A(z - Z2) - A(z - Z3)]} 

X [A,{z - za) A^{z - zi) + A^{z - Z4) A^{z - zi)] 
+ {6 ip2 (2 ifi + ip2) sm[A{z - Z2) - A{z - Z3)] sm[2A{z - z^)] } 
X [A,{z - zi)[As{z - Z2) - A^{z - Z3)] 
xA,{z - zi)[A,{z - Z2) - A,{z - Z3)]] 
+ {2((^i + 2^2)^ -6v?2 (2^1+ v?2) sin^ [A{z - Z4)]} 

x[A;,{z - Z2) - A^{z - Z3)][A^{z - Z2) - A^{z - Z3)] 
+2 (2 <pi + ^2? A,{z - za)A,{z - Z4) (3.18) 

and 

Sij(^--i)(^^j) = 2 I ((/?! + 2 (/72)^ - 6 (/72 (2 </7i + (P2) sm'^[A{z - Z4)] 

-6 sin^(2; - Z2) ip2 (2 </?i + (P2) cos[2 A{z - Z4)] - sin[A{z - Z3)]'^ 

X {ifii + 21P2) - (fi2 (2 ifii + ifi2) sin^[A{z - 24)]) 
+ 3 ip2 {2 (fi + (P2) cos[A{z - Z3)] sin[2 A{z - Z2)] 
X sin[2 A{z - za)]^ Az{z - zi)A2{z - zi) 
+ {2 (2v9i + ^2? - 6 sin^[A{z - ^3)]bi(v^i 
+2 (^2) - </?2 (2</?i + ^2) sm^[A{z - za)]]} 
xA^{z - Z2)A^{z - Z2) 
+ {2 ((^1 + 2(^2)^ - 6(^2 (2^1 + (^2) sm2[A(z - ^4)]} 

xA^{z - Z3)A^{z - Z3) 
+ 2{(2^i + ^2? ]A,{z- za)A,{z - Z4) 
- 6 {ip2{2(pi + (P2) sm[A{z - Z3)] sin[2 A{z - Z4)] } 
X [A,{z - Z2)As{z - Z3) + As{z - Z2)A,{z - Z3)] 
-4{{2^i + ^2fcos[A{z-Z3)] } 

X [A^{z - Z2)A^{z - Zi) + A^{z - Z2)A^{z - z^)] 
+ 6 { sin[A(2; - Z2)] sm[2 A{z - ^3)] 

X [-</7i ((/?! + 2 (/72) + <P2 (2 </?i + </?2] sin^[A(2; - 2;4)]] 
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+ (2v5i + ^2) cos[A{z - Z2)\ sm.[A{z - 2:3)] | 
X sin[2 A{z - z^)] [A,{z - zi)As{z - Z2) + A,{z - zi)A,{z - Z2)] 
z — Z2)] + 12 If 2 (2 fi + 
X sin[y4(2; — z^)] — cos[A{z — z^)] cos[A{z — z^)] 
X sin[A(2; — Z2)] + cos[A{z — Z2)] sin[A{z — Z4)] | 
X [A^{z - zi) Az{z - Z3) + A^{z - zi) A-Xz - 2:3)] 
+ 4{(2(/?i + ^2? sm[A(2; - Z2)\ sin[A(2; - z^)] } 

X \A,{z - Zi)A,{z - zi) + A-,{z - z,)A,{z - z^] (3.19) 



Again it can be checked that the expression for 8i{^j){^^j){^^i) in eqn.( |3.18| ) goes to that 



of a single vortex-antivortex pair of type i in the hmit Z2 ^ z^. However, there does 
not exist any similar simple check on the expression Sij(-i)(^-j), involving limits as some 
of the points zi, .., 2:4 are coincident. This is because the orientation of 5*0(3) elements 
juxtaposes non-commuting elements. For the same reasons, has the most 

complicated form of the three orientations ^[f , ^if- In Fig. 5a, the energy density 
Si(_-i)j(^-j) has been plotted for A = 2, a = — 1 where the vortices are placed at the points 
(-3,3), (3,3), (-3,-3), (3,-3) in the plane. For comparison. Fig. 5b shows a plot of the 
energy density £^j(_j)j(_j) for abelian i-type vortices, at the same positions, interacting 
Coulombically. In both plots we have approximated (fi and 1^2 for 1 2; — 2:2! < fc^Y 
the trial functions given in eqns. (p.l8|) . It is clear that energy is localized at the centre 
of each vortex, which is further evidence that core regions exist. Also apparent in the 
nonabelian case is a suppression of the energy density of j type with respect to i type 
vortices, which is a consequence of the (sine)^ terms in eqn.( |3.15|) . 



In the next section we shall outline a method for computing the the energy of 
these 3 basic configurations, from the corresponding energy density. As we shall see, 
this demonstrates that the interaction energy between non-commuting vortices is quite 
different and rather more complicated than the simple logarithmic Coulomb potential 
that appears in the abelian case. 



4 Interaction potential between i and j-type 
vortices 

In this section we shall present a method for calculating the potential energy of the 
4- vortex configurations, whose energy densities were given in eqns.( pTT^ ), ( p.l8| ) and 



20 



( |3.19| ) of the previous section. This involves deriving a set of differential equations in the 
variables zi, ...,Z4 . For definiteness we shall consider the simplest such configuration 
( |3.14| ), with energy density Sii^^i)j(^^j) . The total energy Ei(^^i)j(^j) = J d'^z Si(-i)j(^j) 



can be written as 

^2, Z3, Z4) = Ell + E22 + E12 + E21 



where 



+ 7T < {ipi + 2if2y>\n r' ^'^ ) + 2Et 

+ TT < (2(^i + (^2)'>ln(^^^) + 2E,^ (4.1) 



Ell = -6 < {2ipi + ip2)> J d^z{A,{z - zi)A-,{z - zi) 

X si\i^{A{z - z^) - A{z - Zi))} 

E22 = -6 <^i (2 ipi + ip2)> J d^z{A^{z - Z2)A^{z - Z2) 

X sin^ {A{z - Z3) - A{z - Z4))} 

E12 = -6 <ipi (2 ipi + ip2)> J d^z{A^{z - zi)A-^{z - Z2) 

X ^m^i^Ai^z - z^) - A{z - za))} 

E21 = -6 <(pi (2 (pi + (p2)> J d^z{A^{z - Z2)Az{z - zi) 

X sin^iAiz- Z3) - A{z- Zi))} (4.2) 



The last four terms in eqn.( |4.1|) represent contributions from isolated vortex-antivortex 
pairs of type i and j respectively. This formula also takes into account contributions 
from within the core regions, which gives rise to the core energies E'f and Ej in this 
equation. En, ...,-^21 on the other hand, represent interactions amongst those pairs. 

Using the relation 922(1/(^1 — -22)) = dz2dziA{zi — Z2) = 5^'^\zi — Z2), one can derive 
the following first order differential equations for E12 and £'21 

_ _ 6 < ^. (2^. + > ^^_^^ _ _ _ 

[Zl - Z2) 

dz,E2i = -^^^i^^^i±^sin^(A(.,-.3)-A(.,-.4)) (4.3) 

[Z2 — Zl) 

Once we have obtained solutions for E12 and E21, those of En and E22 follow by 
taking the appropriate limit zi —>■ Z2. Consider first the solution for E12, which can be 
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obtained by integrating eqn.( [4.3| ), 



dZ2 



X / . sm\A{z2 - z,) - A{Z2 - z,)) + H,^ (4.4) 

J {z^ — Zo) 



{zi - Z2) 

where the function H12 satisfies 9^2-f^i2 = 0. Let I12 denote the integral in eqn. 
It is straightforward to write this integral in the following manner 
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di^ I {z + uji){z + UJ2) ^ 1 (z + cji)(i + c^2) -i 

Z A {z + U)i){z + UJ2) A {z + UJi){z + U)2) 



where ui = zi — z-^ and uo2 = zi — Z4^. Calculating the integrals in eqn. ( [4. 51) , I12 is given 
by 
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1, / \f o , {^2 - Zi){Zi - Z^) {Z2 - Z^){Zi - Zi) 



-ln(2i -22){-2 



4 I (^2 - 23)(^l - ^4) (^2 - ^4)(^1 - 23) 

(^2 — ^3)(^4 — 2:3) , , . (^2 — ^3)(2;4 — 2:3) , , \ t A c\ 

+ 77^ ~\( \Mz2-Zi) + — — -ln^2-24 (4.6) 

4(2:2 - ^4) (^1 - ^3) A[Z2 - Zi)[Zi - Z'i) 

In a similar manner, we can solve for £"21 

E21 = -6 < (Pi{2lPi + ip2) > l2l{zi,Z2,Z3,Z4) + H21 (4.7) 



with dz^H2i = 0. The integral I21 in eqn.( [4.7|) is given by I12 of eqn.( [4.6| ) with zi ^ Z2. 

What remains is for us to determine the function H12. ( H21 is again equal to H12 
under the exchange zi ^ Z2)- To determine H12 we differentiate £^12 once with dz^ , 
and obtain the differential equation 

^ rr 3 < yi(2(^i + (^2) > (Zl - 2:3) (^1 -^4) (fl-^sKfl-^. 

[Z1-Z2) [Zl - Zz)[Zi - Zi) [Zi - Z2,)[Zi - Zi) 

Solving for H12 one finds 

H,2 = 3 < ^1 (2 ^1 + ^2) > I ' ~ 'i n-zi - ^4) 

L [Zi — Z3j[Z2 — Zi) 

(2:1 — 2:3) (^4 — ^3), /- Iw- -^ 

+ 7 77^ -tH^i - Z3) - -\n{zi - Z2) 

[Zi - Zi){Z2 - Z3) 2 
, (^2 - ^4) (2:1 - 2:3) _ ^ (^2 - ^3) (2:1 - 2:3) _ I 

Defining E*"* = E12 + E21 as the interaction energy between the vortex-antivortex pair 
of type i at {zi, Z2) in the presence of the pair of type j at (2:3, 2:4), one finds 
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{Z2 — Z4) 1^2 — 2:41 (2:1 — ^3) \Zi — Z2\. \Z2 — Z4\ 

{Z2 — Z4J \Z2 — Z^l Zi — Z^ Zi — ^2 , ^2 — ^4 ,,, 

^2 — ^3 n /k2 — ^ak , zi — Z4 \zi — Z2\- , |2;2 — 

Z2 — Z4 JT — Z3 Ih Ih 

+ zi^ Z2] (4.10) 



In eqn.( [4.10| ), we have introduced an infrared cutoff R, which can be taken as 
the macroscopic size of the 2-dimensional system. It will be an important check on 
the total configurational energy Ej(_j)j(_j), that one can take the thermodynamic limit 
R ^ 00 and still obtain finite total energy. To verify this, the various logarithmic 
divergences in the interaction energies given in eqn. ( [4.10| ) have to cancel with the self 
energies En and -E22- The latter are obtained from E12 and E21 in the limit zi Z2 
and Z2 —>■ Zi respectively. Of course in taking this limit one should remember that 
there is an effective short-distance cutoff in the value of l^i — 2:2! which is the core size 
Tc- For separations smaller than this, the fields <fi,<f2 are no longer given by their 
vacuum expectation values, but vanish as — 22I — ^ to maintain the finiteness of 
the energy, as discussed in Section 2. 

With this in mind, it is straightforward to compute the following terms proportional 
to Ini? in Ell + -^225 

-Ini? < Lfi {2 Lfi + Lp2) > {4: + — H \- — 

2 Zi — Z3 Zi — Zs Zi — Z4 

Zi — Z^ ^ Z2 — Z4 ^ Z2 — Z4 ^ Z2 — Z3 ^ 2^2 ~ ^3 j ^^-j 
Zi — Z4 Z2 — Z^ Z2 — Z3 Z2 — Z4 Z2 — Z4 

By comparing eqn. ([4.11|) with eqn.( [4.10D , one sees that indeed the Ini? terms can- 
cel, so that one may safely take the thermodynamic limit. This result is the nonabelian 
generalization of the well known property found in abelian vortex systems, that only 
overall neutral configurations, with equal numbers of vortices and antivortices are rel- 
evant (i.e. have finite energy) in the thermodynamic limit. Neutrality in the present 
context simply means with respect to the group Qg, so that for every vortex of one 



23 



type, we must include the antivortex of the same type. This explains the restrictions 
on the allowed total vorticities in the GCPF of eqn. (|3.11|) in the last section. 

It is interesting to see in the formula for the interaction energy eqn.( 4.lO| ) the usual 
Coulomb terms as well as new logarithm-like terms, which differ from the former by 
the prefactors that are rational functions of the points 2:1,. ..,24. One can, by simi- 
lar methods to those described above, obtain explicit expressions for the remaining 
configurations. Whilst the result for _E'j(_j)i(„j) is very similar to that of 
described above, those of and -E'ij(_j)(_j) are particularly complicated and 

lengthy, although they still involve the logarithm-like functions seen above. Details of 
these expressions will be given elsewhere [12]. 



5 Summary and Conclusion 

In this paper we have made some first steps in the investigation of the statistical me- 
chanical properties of nonabelian vortices in two spatial dimensions. We have done this 
with particular reference to vortices characterized by the fundamental group Qg. An 
ansatz was presented that describes isolated vortices, which as we showed in the ap- 
pendix, does indeed describe configurations of minimum energy. To eventually describe 
statistical mechanical properties, it is necessary to discuss multi- vortex configurations, 
which would be present in a reahstic system such as nematic liquid crystals [9], at 
finite temperature. We computed the energy density of all relevant configurations with 
a maximum of four vortices, and gave a method for calculating the total energy from 
this. 

The kinematic arguments discussed in Section three simplified the discussion some- 
what, since they implied that only two out of possibly four species of nonabelian vortices 
would be relevant to the statistical mechanics of the nonabelian gas. Even so, the form 
of the energy and energy density of the many different 4- vortex configurations are very 
much more complicated than the corresponding expressions for abelian vortices . An 
important check on these results was that the energy of these configurations remained 
finite in the thermodynamic limit, in which the size of the system becomes infinite. 

As stated in the introduction, the results presented in this paper should provide 
a starting point in which to investigate amongst other things, nonabelian screening 
mechanisms which generalize those already investigated for abelian vortices [1,6]. A 
step in this direction would be to derive a nonabelian Poisson-Boltzmann like equation 
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for the linearly screened potential between two test vortices, using the form of either 
the energy density or energy as given in Section 3 and Section 4, [12]. There are many 
additional problems that one could investigate in the context of nonabelian vortices and 
K-T hke phase transitions, which have already been studied in the abelian case. For 
example one could consider the effects of putting the system on a surface of non-trivial 
topology e.g. a sphere [16]. In addition it would be interesting to consider further 
modifications of the nonabelian system by including vortex induced Berrys phases. In 
the abelian case, the presence of siich phases was shown to have a dramatic effect on 
the nature of the K-T phase transition [17] . 
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Appendix 



A 
B 



Stability of nonabelian vortices In Section 2 (see eqn.( |2.1(i|) ) and Table 1, we presented 



an ansatz for describing isolated nonabelian vortices described by the fundamental 
group Qg. Specifically these were of the form Ga{0)^diag{r)Q~^{9), where Qa{(^) ,« = 
i, j, k were rotations about the x,y or z axes by angles 6/2, (with the angles 6/2 replaced 
by —6/2 for the corresponding antivortices ). In this appendix we will give a proof that 
these ansatzes are the physically correct ones in that they describe vortices with lowest 
energy^ To accomplish this we shall simply deform our original ansatz by boosting the 
specific choice of 5*0(3) element described above, by an arbitrary rotation depending 
on three angles a, (3, 7 which will be functions of 6 and r in general, and argue that such 
deformations always lead to an increase in the vortex energy. The only constraint that 
the boosts have to satisfy is that they preserve the boundary conditions satisfied by the 
group elements Qa as discussed in Section 2 of the paper (see eqn.( |2.1lD ) This simply 
guarantees that the closed curves in M which Qa parameterize, are homotopically 
invariant under the deformations. For example, a generalized ansatz for vortices of 
type k will be 

where in eqn.( |B.lD a, j3 are the single valued function of 6, while 
^{6 = 0) = 7(6' = 27r) + TT , in order to preserve the boundary conditions as discussed 
above. In showing that perturbations of our original ansatz of eqn.( p.lOD always in- 
crease the energy of a vortex, we shall restrict ourselves to regions outside of the vortex 
core Tc, where the fields (pi{r) and ^2{f) can be approximated by their vacuum expec- 
tation values. Although the particular numerical solutions described in Section 2 used 
our original (unperturbed) ansatz for Qa{6) , the property that they asymptotically go 
to their vacuum values should be true even for a modified ansatz, (see the discussion 
after eqn.( |2.18| ) of Section 2). Outside vortex cores, the energy density £ijk of the field 



configuration taken the form 

^We will do this explicitly for fc-type vortices; the calculation for i and j types follow in a very 
similar manner with the same conclusions concerning stability. We will therefore only present the 
calculation for fc-type vortices only. 
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ijk 



<^i^2 > {2a'^ + 6 cos{2 p) a'" + 8 (3'^ - 8 sin(/?) a'7' - 47 



/2 



3/2 



/2 



< ^2 > 



7a'^ cos{2 13) a" ^ 3 cos(2 {(3 - -f)) a'^ ^ 3 005(27) a 



/2 



2 2 

^ 3cos(2 (/3 + 7))«^' _ 3 g.^^^ _ 2 ^) + 3 sin(/? + 2 7) a' /3 



4 



+5 p'^ - 3 cos(2 7) p'^ + 4 sm(/?) a' 7' + 2 7 



3/2 



/2 



+ 



< ^1 > 



7a'^ cos(2/?) a" 3 cos(2 (/3 - 7)) a'^ 3 005(27) a 



/2 



2 2 

,/2 



3 cos(2 (/3 + 7))a' 



+ 3 sin(/5-27) a'/?' -3 sm(/? + 27) a'/?' 



+5 f3''^ + 3 cos(2 7) f3'^ + 4 sin(/3) a' 7' + 2 7'^ 



(B.2) 



where in eqn.( [B.2| ) '/'indicates differentiation with respect to coordinate 6 and the 
vacuum parameter a. Using the relation between < (pi > and < ip2 > 



< V^2 > = 0, < (fii > 

we rewrite the energy density of eqn. (|B.2| ) as 
where 



Mi. 



( hi qi ri \ 

<?i gi 

\ri Pi/ 



(B.3) 



(B.4) 



(B.5) 



f = a' , , 7' and 



hi = -[(14 + 8a + 14a2 + 6 cos(2/3) + 24a cos(2/3) + 6a2 cos(2/5) 

- 3 cos(2 {/3 - 7)) + 3 cos(2 {f3 - 7)) - 6 cos(2 7) + 6 cos(2 7) 

- 3 cos(2 (/3 + 7)) + 3a' cos(2 (/? + 7))] 
= [5 + 8a + 5a^ + 3 cos(2 7) - 3a^ cos(27)] 

pi = 2 {a -if 

gi = 3 (a^ - 1) cos(/?) sin(2 7) 

ri = 2 (a - 1)^ sin(/5). (B.6) 
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Now, if one can show that all of the energy eigenvalues of this matrix are positive, 
it would be sufficient to show the stability of our original ansatz. The characteristic 
polynomial C{x, a, /5, 7, a) associated with Mj^is given by 

-X^ i9i + hi+Pi)-X [-91 h - giPi- hi pi + gi^ + ri^] 

+Piqi'^ + giri"^ - gihipi = C{x) 

(B.7) 

To analyze the possible signs of the eigenvalues x (which are solutions of the equation 
C{x) = ), we shall first study the extremal points of C{x) with respect to the variable 
X- The general form of C{x) is illustrated in Fig.4. If we can show that the values of x 
at the two extremal points ( denoted by Xei, Xe2 in Fig.4 ) are bounded below by zero, 
then it follows that at least two of the three eigenvalues are also similarly bounded. 
To show that the remaining eigenvalue is > 0, we need only check the value of C at 
X = 0, is < 0. 

The values of x at extremal points of C are given by 

3x^ -2x[gi + h+pi] + {gipi + hpi+gihi-qi^ -n^) = (B.8) 

This has two roots given by 

X = l [{gi + hi+ pi) + ^ {gi + hi+ piY - 3{gih + giPi + hpi -ql- rl)] (B.9) 

To get real and positive solutions for x in ( |B.9|) , requires the following two con- 
straints 

gi+Pi + hi > (B.IO) 
gihi + giPi + hipi - ql - rl > (B.ll) 



As we shall see below, constraint (|B.10| ) is satisfied because each of the functions 



Pi, hi and gi are separately positive. From the definitions of eqn.( |B.6| ) pi is manifestly 
positive, so we will study gi and hi instead. Consider first gi. To find a lower bound 
on this function for arbitrary angles a, (3, 7, we consider varying the vacuum parameter 
a to find out if gi has a minimum value. We find 

^ = 8 + 10a-6a cos(27) (B.12) 
oa 

g 

(^ext = — TTTT (B.13) 

10 — 6 cos(27) 



da? 



10-6cos(27) (B.14) 
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Obviously, > so the extrema of gi found by varying the parameter a is a 

minimum and 

9sin(2 7)' 

5-3cos(2 7)- ^ ^ 

Hence it is clear that gi > 

Moving on to the function hi, we can apply the same ideas and find that 

dh 3 

= 2 + 7a + 6 cos(2/?) + 3a cos(2/5) + -a cos(2 (/3 -7)) 

3 

+ 3acos(2 7) + - a cos(2 (/? + 7)) (B.16) 



-4-12 cos(2/5) 



14 + 6 cos(2 /3) + 3 cos(2 {(3 - 7)) + 6 cos(2 7) + 3 cos(2 {(3 + 7)) 
d'^hi ^ ^ 3 cos(2 (/? - 7)) ^ ,^ , 3 cos(2 (/3 + 7)) 



7 + 3 cos(2 (3) + ^ + 3 cos(2 7) + 

(B.17) 



It is easy to see that > so again aext describes a minimum of hi with 



9 cos(/3)^ (18 - 14 cos(2 /?) - cos(2 /3 - 4 7) - 2 cos(4 7) - cos(2 /5 + 4 7)) 

2 (14 + 6 cos(2 /3) + 3 cos(2 (/3 - 7)) + 6 cos(2 7) + 3 cos(2 (/5 + 7))) 

(B.18) 



hi |„ 



By inspecting the numerator in eqn. (p.l8| ), one can convince oneself that again 
hi |mm> 0. Thus we have shown that constraint ( [B.10| ) is satisfied. 

We use the same ideas as above, to prove that eqn. (|B.ll|) is always satisfied. It will 
be sufficient to prove the bound 

hiPi-ql-rl > (B.19) 

since we have already shown the positivity of hi,gi and pi. In fact, after some manip- 
ulations one can simplify the expression for hipi — q\ — r\ and obtain 

hiPi -ql-rl = 2(a - 1)^ (^^lA^ (_i + « + 3 cos(2 7) + 3 a cos(2 7))^^ (B.20) 

which is now in manifestly positive form. Hence the second constraint eqn. ( p.ll|) is also 
satisfied, and so, as explained earlier, we are guaranteed that at least two of the three 
eigenvalues x ^"^^ positive. To prove that the third eigenvalue is itself non-negative, it 
is sufficient to check the value of the characteristic polynomial C(x) at x = 0. We find 

C(x) 1^=0 = -{gihipi) +piqi^ + giTi^ 

= -8{a-l)^a + 2f{2a+iycos\(3) (B.21) 
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which is a negative semi-definite quantity, so that we have the situation as depicted 
in Fig.4 where all three eigenvalues are positive. This completes the proof that all 
perturbations of the original ansatz describing isolated nonabelian k-type vortices in- 
crease the energy. The vortices described by this ansatz are therefore stable. As stated 
earlier, a similar analysis applied to i and j type vortices shows that the corresponding 
ansatz also describe stable defects. 
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Table Captions 
Table 1 : 

Table 1 lists the SU{2) group elements g of isolated nonabelian vortices corre- 
sponding to elements h, of the discrete fundamental group Qg, which are also 
given. In addition, the associated SO (3) group elements Q are listed. 

Table 2 : 

Table 2 shows how the similarity transformations generated by the matrices U : 
Qi — U Qi , Qj — U Qj , act on the SO{S) group elements of i and j 
type vortices. Also given is the action of U on the fields ^diag- 
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Figure Captions 



Figure 1 : 

Figure 1 sliows tlie ellipse described by the equation Lpf + + <fi <f2 = (a+v) 
for various values of A, A', that corresponds to the vacuum expectation values < 
^diag > of the modcl considered in the text, when p = 0. When p ^ 0, < ^diag > 
takes on discrete values only, given by the intersection of the ellipses with the 
dashed radial lines in Fig.l. At these values of ^diag as well as those given by 
the intersection of the dotted radial lines with the ellipses, there is an unbroken 
U{1) symmetry. 

Figure 2 : 

Figure 2 shows numerical solutions (dotted curves) and trial solutions (solid 
curves) for the fields (fi and (p2 which minimize the energy of an isolated vortex of 
type i, for A = 2, 0.01; A' = and a — —1. Also indicated is the approximate 
position of the core region r < Tc- 

Figure 3 : 

Figure 3 illustrates the relative positions of the regions Ijj, Ijk and I^i of stable 
ij, jk and ki type vortices, with respect to the allowed values of < ^aiag > in the 
case p — O.Boundary values of the vacuum parameters a = < (fi2 > / < (fii > 
are also given. 

Figure 4: 

Figure 4 illustrates the form of the characteristics polynomial C{x) defined in 
appendix A, as a function of x- Xo, Xi and X2 (crosses in Fig.4) are the eigenvalues 
corresponding to the solutions of the equation C(x) = 0. Xei, Xe2 denote the 
values of x at the turning points of C(x). 

Figure 5: 

Figure 5a is a plot of the energy density of four non abelian vortices 

at positions (-3,3), (3,3), (-3,-3) and (3,-3) in the x — y plane, with A = 2, 
A' = and a — -1. For comparison, Fig.5b illustrates a plot of the energy 
density £i corresponding to abelian vortices placed at the same position. 
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/i e Qs ge SU{2) G e 50(3) 



1 

^ = ( n ) 5i = ( " ~ -ul+e, ) ft = I cos(| + 20o) sin(i + 26,) 



1 0\ MO 



i \ _ /e'(^°+3) 

° V e ^ 4V -sin(| + 2^o) cos(|+20o) 

-A _/cos(^o + !) -sin(^o + !)\ /cos(f+2^o) sin(f + 2^o) 

1 Oj Uin(^o + !) cos(^o + !)j l,_3in(i+2.o) cos(|+2.o) 

^ ~^\ / cos(^o + |) ^.sin(^o + |)\ , / cos(f + 2^0) sin(J + 2^o) ^ 



" V 4 y 1 ^ sin(|+20o) cos(f + 20o) 



1\ _i_ / cos(0o + |) sin(^o + f: 



-10; \^_sin(^o + |) cos(0o + f) 



cos(|+26'o) -sin(| +26*0)" 

a/' = I 1 



\sin(f+26lo) cos(| + 26lo) 



•n x / //) , e^ ■ w/3 , / cos(| + 26»o) -sin(| + 26io) 0' 

oj =iisin(^o + !) cos(^o + !)j ^.=[--(2+2^0) cos(f+2^o) 



1 y vol 



-1 \ /e*(^''+ 




-I) 9-r=[ Q ^-^i0oH^] G-i=\0 cos(^ + 2eo) sin(0 + 2eo) 




+ 

- sin(6' + 260) cos{e + 26q) ^ 



-a = 



/-I 0\ ^/cos(^o + f) -sin(^o + f)\ . ^/cos(^ + 2^o) sin(^ + 2^o) 

(,0 -ij Uin(^o + f) cos(^o + f)j Usin(^ + 2^o) cos(^ + 2^o), 

-1 0\ / cos(^o + f) -^sin(^o + f)^ ^ "ttf."! 0^ 

-ij = U-in(^o + f ) cos(^o + |) j ^-1=1— (^^+2^0) cos(^ + 2^o) 



Table 1: 
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Figure 1: 
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